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Abstract 

In data taken with the NA48 experiment at the CERN SPS in 1999, 730 candidates 
of the weak radiative hyperon decay A7 have been found with an estimated 

background of 58 it 8 events. Prom these events the rP A7 decay asymmetry 
has been determined to q;(E!° A7) = —0.78 ib O.lSstat ± 0.06syst, which is the first 
evidence of a decay asymmetry in — > A7. The branching fraction of the decay 
has been measured to be Br(H° A7) = (1.16 ± 0.05,tat ± 0.06,y,t) x lO'^. 
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1 Introduction 



The S° — > A7 decay asymmetry plays an important role in solving a long stand- 
ing discrepancy between the Hara theorem [1] and the observed decay asymmetries 
of weak radiative hyperon decays [2,3]. The Hara theorem states that the parity- 
violating amplitude of weak radiative hyperon decays vanishes in the SU(3) limit. 
Accordingly, the decay asymmetries will vanish in this case. Introducing weak break- 
ing of SU(3) symmetry one expects to observe small decay asymmetries. In contrast 
to this, a large negative decay asymmetry in the weak radiative decay — > jrf 
was first measured by Gershwin et al. [4] and later confirmed [5]. To address this 
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observation, models were developed which tried to obtain large decay asymmetries 
in spite of weak SU(3) breaking. One category consists of pole models, which sat- 
isfy the Hara theorem by construction, and approaches based on chiral perturbation 
theory. They predict negative decay asymmetries for all weak radiative hyperon de- 
cays [6,7,8]. Vector meson dominance models and calculations based on the quark 
model on the other hand violate the Hara theorem. It has been shown that the Hara 
theorem is generally violated in quark model approaches [9]. This second group of 
models favours a positive decay asymmetry for the channel S° — > A7 [10,11]. There- 
fore, the decay S° — > A7 plays a crucial role in differentiating between the groups 
of models. The only previous measurement of the S° — > A7 decay asymmetry has 
been performed at FNAL, giving a value of —0.43 ± 0.44 [14] ^ , and is not able to 
make the distinction. 

In addition to the decay asymmetry, the branching fractions of radiative S° decays 
are similarly difficult to calculate. Predictions on the 'EP — > A7 decay rate vary by 
about one order of magnitude [6,7,8,10,12,13]. Two experiments have reported a 
branching fraction measurement [14,15], with the measured values being about two 
standard deviations apart. The latter measurement is an earlier NA48 measurement 
based on a data set of the year 1997, which contained only 31 S° A7 candidates. 

In this letter we report a new measurement with greatly enhanced precision on both 
the 5° — > A7 decay asymmetry and branching fraction. 

2 Experimental setup 

The NA48 experiment has been designed for the measurement of direct CP viola- 
tion in the decay of neutral kaons. However, due to the comparable life times, decays 
from A and rP hyperons could also be recorded by the NA48 detector. The experi- 
ment featured two almost collinear beams of neutral kaons which were derived from 
450 GeV/c protons hitting two fixed beryllium targets of 40 cm length and 0.2 cm 
diameter, which were located 6 m (near target) and 126 m (far target), respectively, 
upstream of the decay region [16]. The present measurement is based on data from 
a special two-day run period in 1999 where only the near target was in operation. 
The proton beam was extracted from the CERN SPS in a 2.4 s spill every 14.4 s 
with the intensity increased to about 5 x 10^ protons per spill, which is a factor of 
~ 200 higher than in normal running conditions. The near target was located 7.2 cm 
above the primary proton beam axis, which points to the detector centre. It was fol- 
lowed by a sweeping magnet with tungsten-alloy inserts in which the protons not 
interacting in the target were absorbed. The neutral beam was defined by a 0.36 cm 
diameter collimator 480 cm downstream of the target. It had a production angle 
of 4.2 mrad and a downward inclination of 0.6 mrad with respect to the incident 
proton beam, in order to have the axis of the neutral beam pointing to the centre 

^ The original publication [14] quoted a wrong sign of the asymmetry. This has been 
corrected in an erratum in 2002, but not yet in the PDG value [5]. 
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of the electromagnetic calorimeter. The decay region was contained in an evacu- 
ated 89 m long and 2.4 m diameter steel tank, closed by a 0.3% radiation lengths 
thick polyamide (Kevlar) composite window, from the centre of which a beam pipe 
traverses the following main NA48 detector. 

The main NA48 detector elements are a magnetic spectrometer and a liquid-krypton 
electromagnetic calorimeter. The spectrometer is composed of two drift chambers [17] 
upstream and two downstream of a dipole magnet. The magnetic field is directed 
vertically and produces a 265 MeV/c transverse momentum kick. The chambers 
provide 100 /im spatial resolution on the track coordinates, resulting in a total mo- 
mentum resolution oi a{p)/p — 0.6% for 45 GeV/c particles. 

The liquid-krypton calorimeter (LKr) measures the energies, positions, and times of 
electromagnetic showers initiated by photons and electrons. Around 20 t of hquid 
krypton at 121 K are used in a ionization detector. The calorimeter has a structure 
of 13212 square read-out cells of 2 x 2 cm^ cross-section and 127 cm length (27 
radiation lengths) each. The cross-section of the active volume has an octagonal 
shape of 240 cm inscribed diameter. The read-out cells are formed by longitudinally 
stretched copper-beryllium ribbons which act as electrodes. The energy resolution 
of the calorimeter is a{E)/E = 9%/£;e3.2%/v^e0.42% with E in GeV; the time 
resolution is better than 300 ps for photon energies above 20 GeV. 

Further detector elements were used in the trigger decision for the present measure- 
ment: a steel-scintillator sandwich calorimeter (HAG) with a length of 6.7 nuclear 
interaction lengths follows the liquid-krypton calorimeter and measures energies and 
horizontal and vertical positions of hadron showers; a hodoscope composed of 64 hor- 
izontal and 64 vertical scintillator strips is located about 1.6 m upstream of the LKr 
calorimeter. A more detailed description of the experimental set-up can be found 
elsewhere [18]. 

The trigger decision for neutral hyperon decays was based on information from the 
detector elements described above. A positive pre-trigger decision required at least 
one coincidence between a vertical and a horizontal hodoscope scintillator strip, a 
two track signature in the drift chambers, and the sum of the energies deposited in 
the electromagnetic and the hadron calorimeter larger than 35 GeV. The next trigger 
level used information from a preliminary track reconstruction: In order to reject 
Ks 7r+7r~ background events, a mass cut around the kaon mass for the vr+vr" 
hypothesis and a cut on the ratio p+/P- between the positive and negative track 
momenta were applied. This ratio is large for A — > p7r~ decays where the proton 
carries the major fraction of the initial momentum, as opposed to background kaon 
decays where the two charged particles typically have a much lower momentum 
ratio. The trigger rejected a window with 1/5 <p+/p_ < 5 (1/4 and 4 for a brief 
period at the start of the run) and [m^+jr- —mK\ < 30 MeV/c^ (25 MeV/c^), where 
TTiK is the nominal kaon mass. Because of the high rate the triggered events were 
down-scaled with different factors between 4 and 8 during the two-day run period. 
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Fig. 1. Two-track invariant mass with tt'^tt" assumption versus the momentum ratio of 
positive to negative tracks. Only trigger cuts are apphed. The curved structure starting 
at p+/p- — 3 consists of A — > p7r~ decays. 



The plot of mTT+TT- versus p+/p- for the accepted events is shown in Fig. 1. In total, 
approximately 150 million triggers were collected, out of which 18 million events 
contained A hyperons. 

The events were selected by a trigger optimised to select hadronic events for which 
the accuracy of the electromagnetic calorimeter information was not necessary. For 
this trigger, showers from the LKr were read-out only if their seed energy exceeded 
1 GeV instead of the usual 100 MeV. Therefore a cluster energy correction was 
applied, which was determined from the reconstructed invariant 7r° mass in — >■ 
7r"'"7r~7r° events from the same data sample. The correction was energy-dependent 
and ranged between 3% and 0.5% for photon energies between 5 and 20 GeV. 

A Monte Carlo simulation based on GEANT 3.21 [19] for the LKr calorimeter and 
a fast track simulation was used throughout this analysis. Special care was taken to 
adjust beam, target, and collimator geometries and positions to the data. 



3 Method of the Asymmetry Measurement 

For the asymmetry measurement we exploit the well-known decay asymmetry of the 
A p7r~ decay. The A hyperons are longitudinally polarized by the parent process 
S° — s> A7 with a polarization of q;(S° A7) in their rest frame. Effectively, one 
measures the distribution of the angle 0a between the incoming S° (corresponding 
to the outgoing A direction in the S° rest frame) and the outgoing proton in the A 
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Fig. 2. The decay angle 6a in the A rest frame, as measured for both the — A7 and 
the 2° A7r° decay. 

rest frame (see Figure 2): 

, "^^^ = A^o (1 - a(A ^ pn-) ^ A7) cosGa) 
a cos Ba ^ ^ 

In this way the A is polarized by the S° decay and analysed by its own decay into 
p7r~. The minus sign is purely conventional and arises from the fact that the photon 
carries spin 1, which leads to an opposite A spin to that in the process 5° Att^ [20]. 

For cahbration purposes and as a cross-check, we have also analyzed the decay S° — > 
Att", for which the decay asymmetry is well measured. Obviously, when replacing 
A7) with Avr''), there is no difference in the topology for the decays 

S° A7 and S° Avr^ and the definition of the angle Oa is similar. However, as 
explained above, the spin nature of the tt" leads to a sign flip for the longitudinal 
A polarization. 



4 Event Selection 

The A hyperons were identified by two oppositely charged tracks in the spectrome- 
ter. The tracks were required to be at most 4.5 ns apart in time and to have momenta 
of at least 30 GeV/c for the positive and 5 GeV/c for the negative track. For both 
tracks the radial distances to the beam axis at the longitudinal positions of the drift 
chambers had to be greater than 12.6 cm to ensure full efficiency of the chambers. 
The extrapolation of the two tracks to the point of closest approach defined the 
A decay vertex. The distance of closest approach had to be less than 2.2 cm. The 
reconstructed A hyperon was required to have a minimum momentum of 57 GeV/c 
and a reconstructed invariant mass within 2.7 MeV/c^ of the nominal A mass, corre- 
sponding to a window of ±3 standard deviations. Finally, the longitudinal position 
of the A decay vertex had to be at least 6 m downstream of the target, corresponding 
to the end of the final collimator, and at most 40 m from the target. 

Photons were detected by clusters in the electromagnetic calorimeter. The cluster 
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energy had to be above 15 GeV. To ensure negligible energy loss, all clusters had to 
be at least 15 cm and at most 120 cm from the beam axis, well within the sensitive 
volume of the LKr calorimeter, and had to have a distance greater than 2 cm to 
any dead calorimeter cell. To avoid effects from hadron shower contaminations, the 
photon clusters were required to be separated by at least 25 cm from the impact 
point of any track. The cluster-time had to be within ±3 ns of the average of the 
track times. 

The energy centre-of-gravity at the longitudinal position of the calorimeter was 
defined as 

^cog = ^riEi / ^Ei, 

i i 

with the sums running over tracks and photon clusters. The positions are either 
the cluster positions in case of photons or the coordinates of the tracks extrapolated 
from before the spectrometer to the LKr calorimeter. The energies Ei are the photon 
cluster energies or are derived from the track momenta using the nominal particle 
masses. To avoid missing energy in the event, the position of the energy centre-of- 
gravity was required to be less than 7 cm from the detector axis. The target position 
and the energy centre-of-gravity Vcog defined the reconstructed line-of-flight of the H°. 
The distance of the A decay vertex from the S° hne-of-flight had to be less than 2 cm 
to ensure a well- reconstructed event. The decay vertex was reconstructed as the 
point of closest approach between the S° line-of-flight and and the A line-of-flight, 
which was determined from the reconstructed A decay vertex and 3-momentum. 
The longitudinal position of the decay vertex was required to be at least 5 m 
downstream of the target, 1 m before the end of the final colhmator. 

Background from Ks — > 7r"'"7r~7 decays was completely discarded by rejecting all 
events with an invariant mass, under a 7r^7r~7 hypothesis, within 7.5 MeV/c^ of the 
nominal kaon mass (corresponding to 2.5 standard deviations of the resolution). Ad- 
ditional background comes from — > tt'^tt^tt^ decays as well as from time overlap 
of two events or hadronic interactions in the collimator region. Both accidental and 
hadronic background contain a real A decay together with one or two photons from 
a 7r° decay. In order to reject this background together with Kl 7r"'"7r~7r° decays, 
we rejected all events which had one or more additional LKr calorimeter cluster 
within 5 ns of the average track-time and more than 25 cm from the track impact 
points in the calorimeter. Decays of S° — > Att^ with one lost photon are strongly 
suppressed by kinematics to enter the signal region, even if resolution effects are 
taken into account. 

The A7 invariant mass distribution of the candidates is shown in Figure 3. Within 
±7.8 MeV/c^ (3 standard deviations) of the S° mass 730 S° A7 candidates were 
reconstructed. As all remaining background sources have a flat distribution in a close 
invariant A7 mass interval, as seen from Monte Carlo simulation, we have used mass 
side-bands from 1.3 to 1.305 GeV/c^ and from 1.325 to 1.335 GeV/c^ for background 
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Fig. 3. Distribution of the invariant mass for A7 decays. The signal region 

is indicated by the vertical arrows. Also shown are the side-band regions used for the 
background subtraction. The sharp increase of events below 1.3 GeV/c^ is due to H*^ Avr^ 
events with a lost photon. 



subtraction. From this, the background was estimated to be 58.2 ±7.8 events, which 
corresponds to (8.0 ± 1.1)% of the signal sample. It is dominated by hadronic events 
with some remaining contributions from accidental background and mis-identified 

Kl 71^71^71^ decays. 

For the selection of the normalization channel S° —>■ Att^ only a few changes were 
applied with respect to the S*^ —>■ A7 selection. At least two photon clusters within 
5 ns and with energies above 6 GeV were required. The invariant mass of the two 
photons, computed using the cluster energies and positions and the reconstructed 
S° decay vertex, had to be within 3 standard deviations (about 10 MeV/c^) from 
the nominal vr'' mass. No cut on additional LKr calorimeter clusters was applied for 
the S° — > Att^ selection. 

After applying all selection criteria, the total yield of S° Avr" candidates was 
61828 events within a window of ±3.0 MeV/c^ {3a) around the nominal S° mass 
(Fig. 4). Backgrounds are negligible in this channel. 



5 Measurement of the Decay Asymmetry 

As described, the asymmetry is analysed in the angle cos Ba. The resolution on 
cosBa is 0.01, as determined from the Monte Carlo simulation. Fig. 5a shows the 
cos 6a distributions for all S° A7 candidates and for the properly scaled side-band 
events compared to the isotropic Monte Carlo simulation. 

The ratio of data over Monte Carlo corrects for the detector acceptance. To account 
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Fig. 4. Distribution of the invariant 'E^ mass for reconstructed 
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for the background contamination, an effective asymmetry ctbkg of the background 
was determined by fitting tlie side-band events with dN/dcos<d\ oc 1 — atkgCosGA. 
The result of atkg = 0.13 ± 0.32 (together with the background fraction of (8.0 ± 
1.1)%) was used to apply a smooth correction to the cos 6a distribution of the 
signal events for the background contamination. The ratio of background-corrected 
signal over Monte Carlo is shown in Figure 5b. A fit with the normalization and 
the asymmetry a{A p7r^)a(S° A7) as free parameters was performed in the 
range —0.8 < cosGa < 1, where both data and Monte Carlo statistics are large. 
The result of the fit is a{A p7r~)a(S° — > A7) = —0.50 ± 0.12^1^^ with a P^r 
degree of freedom of 9.4/7. After dividing by a(A —>■ pn^) = 0.642 it 0.013 [5] we 
obtain a(S° ^ A7) = -0.78 ± 0.18,taf 

Several contributions to the systematic error have been studied. The largest un- 
certainties come from the effective asymmetry of the background (±0.052) and 
the background normalization (±0.019). Additional uncertainties may arise from 
a possible incorrect modeling of the detector acceptance. This was investigated with 
S*^ —>■ Att^ events, which have sufficient statistics, and which are, due to the lower 
Q value, more sensitive to possible inefficiencies near the beam pipe. Different geo- 
metrical acceptance cuts in the beam pipe regions and on the energy as well as 
uncertainties in the beam profile simulation and the pt distribution of the A hyper- 
ons have been studied, resulting in a variation of ±0.007 on the decay asymmetry. 
The effect of transverse S*^ polarization on the decay asymmetry was studied and 
found to be negligible compared to the other uncertainties. Finally, the uncertainty 
on the used value of a(A pvr~) contributes ±0.016 to the systematics. 

By adding the single components in quadrature, the total systematic uncertainty is 
±0.06, which is small compared to the statistical uncertainty of the data (Tab. 1). 
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Fig. 5. The H*^ — > A7 decay asymmetry: (a) cos 0a distributions of signal candidates 
(crosses), scaled side-band events (shaded), and isotropic Monte Carlo events (dashed), 
(b) Ratio of background-corrected signal candidates over isotropic Monte Carlo simulation. 
The line shows the result of the fit. 



An important check for the validity of the H° —>■ A7 analysis is the measurement 
of the decay asymmetry in S° — >■ Att^, where the data statistics is a factor of 100 
higher. Both channels are two-body decays and the decay S° — Avr^, due to the 
lower Q value, is even more sensitive to effects from the beam-line geometry and 
polarization than S° A7. The decays H° Avr" were recorded with the same 
trigger, analysed with the same definition of cosGa, and used a similar Monte Carlo 
simulation as for the H*^ — > A7 analysis. From our data we measure a combined 
S° A7r° asymmetry of a(A p7r-)a(S° Avr^) = -0.257 ± 0.011,^,^. The 
excellent agreement between this result and the best published measurement of 
a(A p'K-)a{~P A7r°) = -0.260 ± 0.006 [22] vahdates the method and serves as 
a systematic check. 

As a further check, a second analysis, which was performed using the same data but 
different Monte Carlo samples, gave consistent results. 
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A a(SO ^ A7) 



Background asyninictry 


1 n n^o 


Background normalization 


± 0.019 


Detector acceptance 


± 0.007 


a(A ^ pir-) = 0.642 ± 0.013 


± 0.016 


Total systematic uncertainty 


±0.06 


Statistical uncertainty data 


±0.18 



Table 1 

Summary of uncertainties on the — ^ A7 decay asymmetry a(S° — ^ A7). 



6 Measurement of the S° — > A7 Branching Fraction 

Using the central value of our new result on the decay asymmetry, the overall 
S° — > A7 reconstruction efficiency was determined to be 5.4% from Monte Carlo 
simulation. For the normalization channel S'^ — > An^ the reconstruction efficiency 
was determined to be 0.58%. It is significantly lower than for the signal channel, as 
the proton receives less transverse momentum and escapes more often undetected 
through the beam-pipe. Subtracting the background and using Br(5° Att'^) = 
99.5% [5], we obtain Br(S° ^ A7) = (1.16 ± 0.05,^^*) x lO'^. 

The systematic error is dominated by the uncertainty on the decay asymmetry, 
which leads to a ±4.8% uncertainty on the geometrical acceptance and therefore 
on the branching ratio. Other systematic uncertainties are the detector acceptance 
(±1.1% of the measured value), the S° polarization (±1.4%), and the background 
estimation (±1.1%), which have been determined in a similar way to those for the 
decay asymmetry measurement. To account for the unknown S° polarization in our 
data we have varied the polarization in the simulation by ±10%, which is similar in 
magnitude to the measured transverse S° polarization of about —10% for a proton 
beam energy of 800 GeV and a production angle of 4.8 mrad [21]. This variation re- 
sulted in a ±1.4% change of the measured branching fraction. An additional external 
uncertainty arises from the knowledge of the S° — > Att" decay asymmetry (±0.5%), 
since it affects the selection efficiency of the normalization channel. Finally, the 
statistics of the normalization channel and of the signal and normalization Monte 
Carlo simulations contribute ±0.7% to the systematic uncertainty. All other sources 
of systematic errors were found to be neghgible. The total systematic uncertainty 
on Br(S° — > A7) is ±5.3%. A summary of uncertainties on the branching ratio 
measurement is given in Table 2. 

A second independent analysis has been performed also for the branching fraction 
measurement and yielded a consistent result. 
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A Br(SO ^ A7) 




[10-3] 


A7 decay asymmetry 


± 0.056 


S° —>■ Att^ decay asymmetry 


± 0.006 


Background estimation 


± 0.013 


Detector acceptance 


± 0.013 


polarization 


± 0.016 


Statistics of MC and 2° A7r° 


± 0.008 


Total systematic uncertainty 


±0.06 


Statistical uncertainty data 


±0.05 



Table 2 

Summary of uncertainties on Br(E!^ — ^ ^'j)- 



7 Summary and Conclusion 

Prom the data of two days of high intensity Kg data taking of the NA48 experiment 
730 'EP — > A7 candidates have been found with an estimated background of 58.2±7.8 
events. By comparison with an isotropic Monte Carlo simulation, the S° — > A7 decay 
asymmetry has been found to be 

q;(S° ^ A7) = -0.78 ± 0.18,,., ± 0.06.,.,. 

This is the first evidence for a negative decay asymmetry in this channel [23]. As 
has been indicated in the introduction, the H° — > A7 decay asymmetry is a crucial 
input parameter for theoretical approaches to weak radiative hyperon decays. The 
negative sign of the asymmetry clearly prefers models consistent with the Hara 
theorem, while, on the other hand, it cannot be easily explained by quark models 
and models using vector meson dominance. 

In addition we have determined the S° — > A7 branching fraction to 

Br(HO ^ A7) = (1.16 ± 0.05.,,, ± 0.06.,.,) x 10-^, 

which is the most precise measurement of this decay rate so far. It confirms the 
previous FNAL measurement [14], while being about 1.9 standard deviations below 
the previous NA48 measurement [15]. 
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